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The field-dependent magnetization m(H, T) of 1- and 2-dimensional 
classical magnets described by the D-component vector model is cal- 
culated analytically in the whole range of temperature and magnetic 
fields with the help of the 1/D expansion. In the 1-st order in 1/D the 
theory reproduces with a good accuracy the temperature dependence of 
the zero-field susceptibility of antiferromagnets % with the maximum 
at T< | Jo 1/D (Jo is the Fourier component of the exchange interac- 
tion) and describes for the first time the singular behavior of x(H, T) 
at small temperatures and magnetic fields: limT^o liniff _^o x(H, T) = 
1/(2|J |)(1 - 1/D) and limjj-o lim T ^o x(H, T) = l/(2|Jo|). 



KEY WORDS: Low-dimensional magnets; magnetic susceptibility; 
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1. Introduction 

A great variety of low-dimensional magnetic systems were synthesized and ex- 
perimentally investigated during the last decades (see, e.g., refs. |l]] and |^]). 
The idealized 1- and 2-dimensional models (without intcrplane or interchain 
coupling and anisotropy) are characterized by a strong short-range order in the 
low-temperature region, whereas the long-range order is ruled out being smeared 
off by the longwavelength spin waves. Complementary to the high-temperature 
series expansions (HTSE, see, e.g., refs. [|[ ||, P), such approaches as the "mod- 
ified spin wave theory" || and the "Schwinger boson mean field theory" § 
were applied to low-dimensional antiferromagnets at low temperatures. These 
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two methods giving very similar results (with a wrong factor in the antiferro- 
magnetic susceptibility in ref. [Q) are not rigorous expansions in the parameter 
T/\ J\ << 1 (J is the exchange integral) but rather some variational approaches. 
The results break down, however, at T ~ | J| and thus cannot describe the sit- 
uation in the whole temperature range. Since the absence of magnetization 
was introduced in ref. Q as some additional self-consistency requirement, the 
generalization for the case with the external magnetic field H is a problem. 

In ref. Q an analytical method of calculation of the physical quantities of 
classical low-dimensional magnets in the whole temperature range was proposed, 
which is based on the 1/D expansion for the model of Z)-component classical 



with the help of the diagram technique developed in ref. (lC[|. For the Heisen- 
berg model (D = 3) in the 1-st order in 1/D the calculated temperature de- 
pendences of the antiferromagnetic susceptibility and internal energy at H = 
turn out to be very good, which is shown by the comparison with the MC sim- 
ulation data for the internal energy of the square lattice (s.l.) classical 
ferromagnet and with the exact results |l2| for a "toy" example of the classi- 
cal linear chain (I.e.) model. In particular, for both models the characteristic 
maximum of the antiferromagnetic susceptibility at T <\Jq\/ D is in contrast 
to refs. § and [| well reproduced. The reason for the efficiency of the 1/1? 
expansion even for D — 3 is that it yields the exact results for the thermo- 
dynamic quantities at T — > and reproduces several leading terms of their 
HTSE expansion (see ref. [|]) interpolating thus between these limits in the 
whole temperature range. The applicability of the approach to classical low- 
dimensional magnets proposed in ref. ^ is not restricted to the case H = 0, 
and it can be applied to the problem of the singular behavior of the antifer- 
romagnetic susceptibility \ AF (H,T) at low temperatures and magnetic fields, 
i.e. the inpermutability of its limits lirny^o limT-*o an d limT^o lim# , which 
could not be up to now described by other analytical methods. The physical 
reason for such a singular behavior is the following. With lowering temperature 
the system becomes locally ordered, and D — 1 susceptibilities transverse with 
respect to the local sublattice orientation tend to the value l/(2|Jo|) (Jo = zJ, 
z is the number of nearest neighbors), whereas the longitudinal one tends to 
zero. At H = there is no preference direction, and the susceptibility of the 
sample is given by the average over the local sublattice orientations, which re- 
sults in lim T ^ liniff^o X AF \H,T) = 1/(2|J Q |)(1 - 1/D). For the Heisenberg 
model the D-dependent factor makes up the well-known number 2/3. To the 
contrary, for the arbitrary small H ^ at sufficiently low T the lowest-energy 
state with the sublattice magnetizations driven perpendicular to the field H 
and tilted in the direction of H is realised. In this state the susceptibility takes 
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up its transverse value: lirrifl^o limT^o X AF (H > T) = l/(2|Jo|). A quantitative 
description of this effect and the calculation of the magnetization m(H, T) of 
low-dimensional classical antiferromagnets in the whole range of temperatures 
and magnetic fields with the help of the 1 / D expansion is the purpose of this 
work. 

The following part of the article is organized as follows. In Section 2 an 
improved version of the diagrammatic l/D expansion for classical spin systems 
with the magnetic field is developed and the results for the magnetization and 
spin-spin correlation function in the 1-st order in l/D are obtained. In Section 

3 the general l/D-results, which are double integrals over the Brillouin zone, 
are calculated analytically and analyzed for the classical linear chain model, for 
which there is no exact solution in the case with the magnetic field. In Section 

4 the results are converted into the form convenient for numerical calculations 
and the analysis at low temperatures for 2-dimensional systems, and the tem- 
perature and magnetic field dependences of the antiferromagnetic susceptibility 
are represented. In Section 5 some important features of the l/D expansion 
and its applicability to the systems with concrete values of D are discussed. 



2. The l/D expansion 

The physical quantities of ferro- and antiferromagnets described by the hamil- 



tonian (1.1) can be expanded in powers of l/D with the help of the diagram 
technique for classical spin systems ||, ^0). Here the consideration of ref. || is 
improved and generalized to the case H ^ 0. We choose the z axis along the 
magnetic field H, the other (transverse) components of the spin vector m are 
designated by the index a = 1,2, .., D — l. The wavevector-dependent transverse 
susceptibility x± (k) = Xa (k) of a classical spin system is related to the Fourier- 
transform of the spin-spin correlation function S aa (r — r') =<m a (r)m Q (r')> 
via the formula Xx(k) = (3S aa (k.), (3 = 1/T. With the help of the diagram 
technique of ref. g] x± (k) can be represented as 

= AW*) (2.1) 
l-A QQ (k)/?J k 

where A QQ (k) is the compact (irreducible) part of S aa (k) given by the diagrams, 
which cannot be cut by the one interaction line /3Jk- For the isotrope systems 
considered here it is not necessary to write down the diagrams for the magne- 
tization m =<m z >, because m(H) can be determined from ( |2.l| ). Indeed, in a 
transverse magnetic field H± « H the magnetization m is simply rotated on 
the angle 9 = H±/H << 1, which results in the important relation 

X±=X±(Q)=m/H. (2.2) 
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The longitudinal susceptibility can be determined now by the formula 

Xz = dm/dH = X x + H(d X ±/dH) (2.3) 
which is much easier than the direct diagrammatic calculation of x z (k). 



The compact part A QCt (k) in (2.1) can be represented in the 1-st order in 
1/D by the diagram set from ref . [ ^completed by the additional diagrams for 
H ^ 0, which can be estimated and selected according to the same rules. The 
general principle here is that the diagrams with multiple irreducible integrations 
over wavevectors (i.e. those that do not separate into products of independent 
simpler integrals) are small as the corresponding powers of 1/D. Thus, in each 
order in 1/D the complexity of diagrams to be taken into account is restricted: 
in the zeroth order in 1/D (the spherical model) only diagrams with the one- 
loop integration over the Brillouin zone survive, and in the 1-st order in 1/D the 
double integrals over the Brillouin zone appear. The large number of diagrams 
in the case H ^ necessitates, however, an improvement of the method, which 
consists in taking into account some diagrams implicitly with the subsequent 
solution of the corresponding equation for A QCt (k). All the diagrams, which 
contribute to A aQ (k) in the 1-st order in 1/D, are represented in figs. 1,2. Note 
that the renormalizcd transverse interaction lines in fig.l contain the unknown 
quantity A QQ (k) itself, which means implicitly accounting for the additional 
class of diagrams of the type 1 and 2 in the fig. 3 of ref. ||. At H = from all 
the diagrams in fig. 2, a survive only the diagrams 1, 2 and 5, and the last term in 
the Dyson equation for the longitudinal interaction line fig.2,b disappears. The 
wavevector dependence of A QCt (k) is due to the diagrams 5-8 in fig. 2, a. There 
is one more diagram 7' that is analogous to 7 and is not represented in fig. 2, a 
to save the place. Taking into account only the diagrams in fig.l results in the 
self-consistent Gaussian approximation (SCGA), which describes rather good 
the thermodynamics of 3-dimensional ferromagnets ]To| . The analytical form of 
A aa (k) in fig.l reads 

A^k^A^+A^^k) (2.4) 

where AaQ D '(k) is the sum of the diagrams represented in fig. 2, a vanishing in 
the limit D — > oo (see Appendix), and A aa is the renormalized 2-spin cumulant 
average given by ||, [bj 

Aaa = * /9 / d D -\ exp(-r 2 )A QQ (C). (2.5) 



7P 

Here A aa is one of the cumulant spin averages of a general type M 

A (0 - (2 6) 

obtained through the generating function A(£) = lnZo(£), 

Z (0 = const •(-W 2 - 1 )V 1 (e) (2.7) 
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+ (1/D diagrams) 




A/WV = 




Fig. 1: (a) Diagrams for the compact part A aa (k) (see also fig. 2, a); (b) block sum- 
mation of transverse loops for the renormalized cumulant one-site 2-spin average A aa ; 
(c) Dyson equation for the renormalized transverse interaction. 




A/WV 



AAA/V -A/WV 
A/WV = AAA/V 




Fig. 2: (a) Additional 1/D-diagrams for A aa (k); (b) Dyson equation for the renor- 
malized longitudinal interaction; (c) ladder equation for the four-spin correlation line. 
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is the partition function of the D-component classical spin, is the modified 
Bcssel function, 

C = /3(H + mJo) + 2^ 2 r, (2.8) 



r is the (D — 1) -component vector perpendicular to H. The last term in (2 
describes the Gaussian transverse fluctuations of the molecular field H + m Jo 
with the dispersion proportional to l a , which leads to the renormalization of the 



cumulant spin averages described by (2.5) for A aa and by analogous formulas 



for the other cumulants entering Aaa D ^(k) (dashed ovals in figs. 1,2). This 
renormalization results from the block summation of the all one-loop diagrams 
with the transverse interaction in fig.l,b; the quantity l a is given by the integral 
over the Brillouin zone 

l a = lj** ^ (2.9) 

2 7(2^1-A QQ (q)/3J q 1 ' 

where is the unit cell volume (vq = Og for the I.e. (d = 1) a nd s.l. (d — 2), ao 
is the atomic space). The (D — l)-dimensional integral in ( [2.5| ) can be simplified 
taking advantage of the symmetry with respect to the transverse variables and 
using the explicite form 



-V,..(€) = ^(l-p)+B'(0^, <->■"<>' 



where 

B(0 = dA(Z)/d£ = Iz,/2(0/Io/2-i(0 (2-11) 
is the generalized Langevin function and B'(£) = dB/d£. Making in ( 2.1C| ) the 



substitution £ 2 = £^/(D — 1) with £ r = 2l]J 2 r and in ( |2.5| ) the partial integration 
to get rid of B', one obtains 

oo 

l --W^m)l drrDcM ~ r2) ^f CHCI (2 ' 12) 



The formulas (2.4), ( 2.12j ) and (2.9) yield the integral equation for the compact 



part of the spin-spin correlation function A QCe (k) entering the basic expression 

By the expansion in powers of X/D the quantities A aa and Aaa fl '(k) in 
(^J) give rise to the terms starting from the zero- and the 1-st orders in 1/D, 
correspondingly, whereas the expansion of all other diagrams neglected here 
starts from 1/D 2 . Before proceeding with the calculations we make a reference 
to the simplest approach — the mean field approximation (MFA) — in which 
no diagrams with the integration over wavevectors are taken into account. In 
this case l a ^ and C = £ = P(H + mJ Q ), and in (EO) A^ D) (k) => 0, 



G 



A, 



£(£)/£• Now witn tne use of (P3) and (!•!) one S ets the Curie- 



Weiss equation m = -B(£) for the magnetization m, which yields the phase 
transition temperature Tq 1FA — \Jq\/D. The latter has no physical significance 
for 1- and 2-dimensional magnets but can be used as a temperature scale. It is 

iperature 9 = T/T^ 1 
Jo\X- Then the formulas (| 



convenient to introduce the dimensionless temperature 9 = T/T^ IFA , mag netic 

T) and @ 



field h = H/\ Jo | and susceptibility % 
can be rewritten as 



Xx. (k) = 



1 - i/G k A k 



~ l a 1 f dq X c 

l a = ~R = 7777 Vo 



D 29 



27T) 



1 - ^G q A q 



(2.13) 



where Gk = (D /9)A aa ( k), v = ±1 for ferro- and antiferromagnets and Ak = 
Jk/Jo- In the integral ( [2.12] ) the product r £> exp(— r 2 ) is at large D sharply 
peaked at r = ro = (D/2 ) 1 / 2 , whereas changes slowly with r. Using the 

expansion of B(£) ( 2.11 ) for D ^> 1 || one can write 



g(0 

c 



1 



JJ 1 + x- 



g( x ) = 



i 



(2.14) 



where 



-2 



x = 2C,/D = x{f) = U(h + vmY/6 A + 81 a r 



1/2 



i + (i + x 2 y/ 2 

; f = r/r (2.15) 



and evaluate ( 2.12 ) by the pass method. In the 1-st order in 1/D for G 
(D/9)A aa one gets 



G 



+ 



1 



9 1 + [l + 4(/i + t/m) 2 /(9 2 + 8y i /2 D 



A (G) 



where 



A<°> = - 



Idg ld 2 g\ 
4dr + Adf 2 )^ 



(2.16) 



(2.17) 



is the 1/D correction to the Gaussian integral ( p. 12 ) and the derivatives of g 
are calculated with the use of (2.14) and (2.15). The 1-st term of (2.16) also 
contains the 1/D corrections due to the corresponding corrections to m and l a . 

Before further proceeding with the 1/D expansion we consider at first the 
limiting case D — > oo corresponding to the spherical model |13[ . In this case 
the quantity A QQ (k) in ( [2~4| ) can be neglected, and one can re place in ( 2.13 ) 
Gk, G q => G, where G = G(mo, f a o) is given by the 1-st term of ( 2.16| ) with the 
zero-order quantities mo and 



laO 



P(G)-1 . 
20G ' 



P(G) 



f/q 



(27r) d 1 - GA q 



(2.18) 



Here for the square lattice with the n.n. interaction the lattice integral P{G) is 
given by P(G) = (2/n)K.(k) with k = G, K(fc) being the elliptic integral of the 
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1-st kind, and for the linear chain P(G) = 1/(1 — G 2 ) 1//2 . For bipartite lattices 
considered here the integral P(G) is the same for ferro- and antiferromagnets 
and independent of the sign of G. For this reason the sign-factor v is dropped 
in the definition of P(G) ( p. 18 ). Note that for 1- and 2-dimensional systems 
P(G) diverges for G — ► 1, which is the rea son for the absence of the long-range 
order. Eliminating now l a0 and using (2.2) in the form 



too 
h 



G 



1 - vG 



one comes to the equation of state of the spherical model: 

BGP(G) = 1 - ml; 



(2.19) 



(2.20) 



which for h ^ should be solved together with (2.19) (in the general case nu- 
merically). For low-dimensional antiferromagnets (y — — 1) at low temperatures 
(6 <C 1) in the field region where m\ < 1 the equation ( 2.20| ) requires G = 1 and 
P(G) ^> 1. For the square-lattice model for P and its derivative this implies 



P(G) =-hr 

7T \ 1 



G 



P\G) 



1 1 



1 -G 



I -ml 



1 

— 7T~ ex P 



8?r 



tt(1 -ml) 



(2.21) 



i.e. the deviation of G from unity is exponentially small: 
G=l-8ex p r 



(2.22) 



V[2(l-m§) 



2\2l 



For the linear chain model the corresponding result reads 1— G 
Now with the help of (2.19) one gets for the magnetization mo = h/2 with only 
exponentially small corrections in the 2-dimensional case due to ( |2.22j ). The 
latter is valid up to the magnetization saturation point h = 2 (i.e. H — 2| Jo|), 
which corresponds to the spi n-flip field o f 3-d imcnsional antiferromagnets. For 
the fields h > 2 according to ( |2.19|) and (|2.20[) to = 1 and G < 1 : 



m S 1 - (6/2)GP{G); G = l/(h - 1) 



(2.23) 



In the zero field case the antiferromagnetic susceptibility x = X± = Xz = 
G/(l + G) monotonously decreases with rising temperature from the value 1/2 
at T = to at T — * oo, i.e. the spherical model does not describe the maximum 
of the antiferromagnetic susceptibility at 9 < 1. 

Now, returning to the 1/D expansion, one can express the 1/D -correction 
term A^ G ^ in (2.16) through the variables of the spherical approximation: 



A< G > - 2G 



y-i 



p - 1 



\_2y-l 2(2y-l) 2 



(-P-l) 2 (3y~l) 
(2y-l)3 



1 



(2.24) 
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and represent the unknown quantities to, l a and Gk (see (2.13)) in the form 

to = mo + mi/D 

la — laO + lal/D (2.25) 

Gk = G + AGJD 



Here the corrections mi and l a \ can be expressed through AGk with the use of 
( |2 . 1 3| ) and the relation m/h = x_i_(0) = Go/(l — vG ), which results in 

mi AGo 
~~h = (1 - vGf 

and 

1 f dq AqAG q 



(2.26) 



'oil 



26» 7 (27r) d (1 - vGXq) 2 



(2.27) 



Expanding now the 1-st term of (2.16) up to the 1-st order in mi and l a x, one 
comes to the 1/D part of the equation (2.4) in the dimensionless form 



m 2m|y G 2 
AGkH — - 7 AG + - -v 



uh 2y-l 



dq 



A q AG q 



2y - 1 "7 (27r) d (1 - ^GA q ) 2 



A (G) +A a/^) ; (2 . 28 ) 



where the quantity A^ D ^ is the nonvanishing in the limit D — > 00 part of 
13 • (-D/^A^a^^k) (see Appendix). The solution of the integral equation ( 2.28 ) 
has the form 

AG k - AG + M k 

where 



M k 



(i/o) 



-A 



(i/o) 



(2.29) 
(2.30) 



and AGq is given by 



AGo 



(2y-l)(A( G '+A 1/I?) )-G 2 



rfq 



A q Mq 



(2n) d (1 - ^GA q ) 2 



(2.31) 



GP'+P{G)- 



2mly 
1 - vG 

(i/o) 



where P' = dP(G)/dG. Now, calculating the quantity A^ ' ; (see Appendix) 
and introducing the function 



dp 



1 



j9p9p-q.i 



9p 



1 - z/GA c 



(2.32) 



one arrives after numerous cancellations at the final results for AGq and M^: 



AGo = 2G jl - [Gpi + P+2mty)] «d/^g 



TO J/ V 



dq g\ 



G 

2 Vo 



dq 7-q 

(2tt)' 1 f q 



GP'+P- 



(27T) 

1 - j/G 



(2.33) 
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and 



M k = 2G v 



dc l gq ~ gq-k 



(2.34) 



where f q = r q + 2mf j yg <i and r q = dr^/dG. The similar results obtained earlier 
[ fli] , p"5f by another method for the particular ferromagnetic case were used for 
the investigation of the phase transition in 3-dimensional ferromagnets. It is 
interesting to note that the function r q ( 2.32| ) is (like P(G)) identical for ferro- 
and antiferromagnets and has a singularity at G — > 1 and q — > 0. In contrast, 
the quantity f q entering (2.33) and (2.34) has for antiferromagnets one more 
singularity at G — ► 1 and q — > b (b is the inverse lattice vector) due to c/ q , 
which disappears, however, in zero magnetic field (mo = 0). This is a formal 
mechanism responsible for the singular behavior of the susceptibility x(^ T) 
of the low-dimensional antiferromagnets in the limit H, T — > discussed in the 
Introduction. 

Before proceeding with the application of the results obtained to the con- 
crete systems it is worth to note some general properties of the k-dependent 
spin-spin correlation function (see (2.1) and (2.13)) that can essentially simplify 
the consideration in the low-temperature range. In partic ular, for 2-dimensional 
ferromagnets in the spherical limit the quantity G ( 2.22 ) is exponentially close 
to unity at h = and low temperatures, which implies exponentially small gap 
in the spin-wave spectrum. Since this property canno t be changed with taking 
into account 1/D corrections, the quantity AGo in ( 2.29 ) should be also ex- 
ponentially small. This is physically clear and can be confirmed by the direct 
analysis || of th e resu lts obtained. On the other hand, the k-dependent con- 
tribution Mk in (2.29) has not to be exponentially small at low temperatures 
and k ^ 0. In fact, the value Mo determines the 1/D correction to the anti- 
ferromagnetic susceptibility in zero magnetic field |8|, which can be expanded 
in powers of 9 <C 1. Thus, by calculation of such quantities of 2-dimensional 
magnetic systems at low temperatures, which are not exponentially small, one 
can use only the quantity Mk (2.34), being much simpler than the expression 
for AGo (2.33). This means that only k-dependent diagrams for the compact 
part of the spin-spin correlation function A aa (k) should be taken into account 
in the low-temperature range, which is a clear advantage of the diagrammatic 
1/D expansion in comparison with the earlier version Jl4| , 15 . The consider- 



ations above can be extended also on 2-dimensional antiferromagnets at low 
temperatures in the field region h < 2 (H < 2| J |), where the magnon gap is 
exponentially small. Here the quantity Gb in ( 2.13| ) should be exponentially 
close to unity in all orders in 1/D. Consequently, the quantity AGo contribut- 
ing to the magnetiz ation and susceptibility of an antiferromagnet (see ( 2.26| )) is 
given according to ( 2.29 ) by AGo = AGb — Mb = —Mb. In the next sections we 
apply the results of the 1-st order in 1/D obtained above to the analysis of the 
equation of state m{H,T) of 1- and 2-dimensional classical antiferromagnets. 
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3. The linear chain classical spin model 



For the linear chain model = Jk/Jo — cos(aofc), and the integrals (2.32) and 
( 2.3S ) can be calculated analytically. One gets r q = 2P(G)/(2 - G 2 - G 2 \ q ) 
with P(G) = 1/(1 - G 2 ) 1/2 and 



AG 



where 



2G(l-m 2 )(l-G 2 ) 



2^,3/2 



2ra\vG 



5 + 3m 2 , ~ 2mlG 2 
2[l + m 2 (l-G 2 )] 



F (3.1) 



1 - F 



m§(l 



vG) 



1-G 2 



f (1 - vG) - 3(1 - ml) - 2vGml 



F = 



l + m 2 (l-G 2 ) 



TH ( 3 - 2 ) 

[(1 + m 2 ) 2 (l - G 2 ) + 2G 2 m 2 (l - m 2 )(l - vG)] 1 ' 2 

The magnetization mo and the par amete r G of the s pherical model in ( |3.l|) and 
( [3^2] ) are given by the solution of ( |2.19 ) and ( 2.2C ). It can be shown that in 
zero magnetic field the results for the susceptibility x are equivalent to those 
obtained by the expansion of the exact solution |Q up to the 1-st order in 
1/D. In the low field and temperature limit h, 9 <C 1 one has mo = h/2 <C 1, 
G = 1 — 9 2 and hence F = (h 2 + 9 2 )~ x l 2 » 1. Taking into acount the leading 
contribution into AGo given by the 1-st term in square brackets, one gets with 
the use of ( 2.25| ) and ( 2.26| ) the following result 



^ = J = 2 



i + D{-W 



.02)1/2 



0{B 2 ) 



(3.3) 



It can be seen that for h = the susceptibility %j_ decreases with lowering tem- 
perature due to the term 9 in (O) and attains the value x± = (1/2) (1 — 1/D) 
at 9 = 0. If /i 7^ 0, then at = /1 2 / 3 the value of \± attains a m inimum and 
then rises to 1/2 at 6 = 0. Note that the singular term in (3^3) becomes of 
order unity at 9 ~ h <C /1 2 / 3 , which is one more characteristic temperature. 
Such a qualitative behavior of the susceptibility of a low-dimensional classical 
antifcrromagnet is in accord with the physical considerations made in the In- 



troduction. The longitudinal susceptibility \ z calculated with the help of (2.3) 
and (3.3) has the form 



dm 
~dh 



1 



1 

D 



2)3/; 



+ 9 + Q{9 2 ) 



(3.4) 



This expression has the minimum at 9 = 3 1 / 3 /? 2 / 3 ^> h, and also the maximum 
at 9 = 3 _1 / 2 /i 3 / 2 <C h (the 3-rd characteristic temperature) where \z — 1/2 + 
(2/D)(/i/3) 3/2 > 1/2. The susceptibilities x±{h, 9) and x z (h, 9) are represented 
as functions of temperature for some field values in figs. 3, 4. 

An interesting feature of the susceptibility \ z manifests itself in the 1/D 
approximation in the low temperature limit (9 <C 1) in the vicinity of the 
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Fig. 3: Temperature dependences of the transverse susceptibility x± = m/h of the 
I.e. Heisenberg antiferromagnet for different magnetic fields in the 1-st order in 1/D. 
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magnetization saturation point h = 2 (i.e. H = H c — 2\Jq\). In the spherical 
limit, adopting mo = 1 — 5mo and G = 1 — 5G with Sitiq,5G <C 1, one can 



simplify the equations (2.19) and (2.20) for the linear chain antiferromagnet to 
8 1 '{25G) 1 / 2 — 2(5too and 28m^ — 6G = 2 — h, which results in the following 
equation for SrriQ in the scaled form 



x - l/(16x 2 ) 



5m /6 2 / 3 ; x = {2-h)/{26 2 ' 3 ) 



(3.5) 



This equation describes the temperature-induced rounding of the transition 
between the dependences mo = h/2 and mo — 1 in the small field interval 
|2 — h\ ~ 9 2 / 3 . N ow, the X/D correction mi determined for |2 — h\, 6 <C 1 from 
2l) and (O), fO) has the form 



mi 



5m (3 - 3F- 1/2 + y 1/2 )/y ; ^ = I + Sm /SG = 1 + 8x 3 (3.6) 



where x is the solution of (3.5). In the limiting cases one gets from (3.7) and 
([3.6]) for the magnetization m = m + nii/D the following results 



rn = < 



r ^ 

2 

1 - 
1 - 



D2(2-/i) 1 /2 [2(2 -h)f 



2/3 



1 

U V 3 



1/2 



(6 1 / 2 -!) 



[&(h~2)} 1 / 2 



d 2/3 ^ 2 - ft < 1 
?i = 2 

6» 2 / 3 < ft - 2 < 1 



(3.7) 



It can be seen from (^J) that in the field region below the saturation point 
h = 2 the temperature-dependent correction to m is positive. Accordingly, the 
susceptibility \z — dm/dh is greater than 1/2 in this region, but the effect is 
not great. With the use of ( [3.6[ ) one can show that for D = 3 Xz,max = 0.518 
at 2 — h = 2.80 • 6 2 ^ 3 . The field dependences of the normalized susceptibility 
Xz of the 1-dimensional classical antiferromagnet are represented for different 
temperatures in fig. 5. It is interesting to note that a qualitatively similar field 
dependence of the susceptibility with a logarithmic singularity at small fields was 
found in ref. fT(| for the quantum linear chain Heisenberg antiferromagnet with 
S = 1/2 at T = 0. There are no physical comments to this effect in ref. ]lq] , but 
it seems now rather clear that the origin of this low- field singularity of a quantum 
antiferromagnet is also the orientation of sublattices perpendicular to the field, 
the quantum effects playing here the role of some "residual temperature" . 



4. The square- lattice classical antiferromagnet 



For 2-dimensional lattice the integrals in AGq (2.33) cannot be calculated ana- 



lytically. For the convenience of the analysis at low temperatures and numerical 
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Fig. 5: Field dependences of the longitudinal susceptibility of the I.e. classical Heisen- 
berg antiferromagnet for different temperatures. 



calculations we introduce instead of the strongly singular r q ( 2.32 ) the weak- 
singular function H 



V'q = 4 ^0 



dp 



A q GA p A p _ q = 1_ [2p{G) 1 (i _ 



G 7(27r) d (1 -GA P )(1 - GA P _ 



G 2 



(4.1) 

in which the divergences of the integrand at p = and p = q at G = 1 are par- 
tially compensated by the nullification of the numerator. In the longwavelength 
region the function ?/> q has the form ||: 



■ In ■ 



1 



7T X ' 



and its derivative ^ q = dipq/dG is given by 

2 1 



7T 1-G' 

2 , a; 

-^ ln T^G ; 



x = (a <7) 2 < 1 - G 
1 - G « i « 1 

a; < 1 -G 
1 - G « i « 1 



(4.2) 



(4.3) 
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At low temperatures at the corners of the Brillouin zone tpq 
In terms of ipq the function AGo ( 2.33| ) can be written as 



-1 + 0(1-G). 



AGo = 2G{ (l-i/) 



GP' + P+3m 2 y f dq 1 



2y-l 



■ V 



GA 



G 



■ v 



dq a^q 1 + GA q 3y - 1 
(2^Y 1 - aV> q 1 + G q A q ~ 2j/- 1 

;+P(G) 



aV»q 1 

■ V(, 



GqAq 

dq 



('0 



(27r) d 1 - aipq 
dq 1 



1 



GqAq 



1 + GA C 



2 7(2^) d l + GqAq 



{2n) d 1 - o^q 1 -f 
GP , +P(G)-' 2m2 ° V 



GqAq 



(4.4) 



1 - z/G 



where y = 1/(6G), a = G 2 /(2y - 1) and 



Gri — G 



2(1 - i/)mgy 



(2y-l)(l-o^ q ) 
For the quantity Mk ( ^.34 ) one gets in a similar way 



M k - 



2G 



dq (l-G 2 A q )(g q - gq _ k 
2y - 1 "7 (27r) rf (1 - o^ q )(l + G q A q ) 



(4.5) 



(4.6) 



Putting k = b in (4.6) in the antiferromagnetic case iv — —1) and taking into 
account only the exponentially great terms with P'(G) in AGo (4.4) at low 
temperatures in the field range h < 2, one arrives at the result 

dq 1 A q 



AG = ~M b = 4a v 



(27T) d 1 - O^q 1 + GqAq 



(4.7) 



which confirms the conjecture made at the end of Section 2. Since at 9 <C 1 
the quantity a = (6/2) /(I - 9/2) <C 1 and due to the functions Vq are 

integrable ones, one can expand (4.7) in powers of aipq and then of 9 to get 
the development of the 1/D correction to the magnetization mi (2.26) at low 
temperatures. In the lowest order in one gets 



m i = TT 



h9 1 - P(G) 
2~ 



1 



2ml 



(4.8) 



where G = G(l — 2mj}) and too — ft/2. This correction is negative for h < 2 1 / 2 
and positive for 2 1 / 2 < h < 2. In the case h ^ the quantity 1 — G can be inter- 
preted as proportional to the field-induced gap of the out-of-plane spin waves, 
which makes the lattice integral P(G) in (|4.8| ) not divergent at low tempera- 
tures. The more detailed physical interpretation of the spin wave dynamics in 
low-dimensional magnets requires, however, the dynamical generalization of the 
diagram technique used here. In the small field region, where G = 1, with the 
help of ( 2.25D and (2.2) in the 1-st order in 1/D one gets 



m 
~h 



1 

D 



In 



l-G + ft 2 /2 



+ 9 + Q{9 2 ) 



(4.9) 
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where G is given by ( 2.22) ) . In the case h = the ln-term in (4.9) is identically 
equal to —1, and — ► (1/2) (1 — 1/D) in the limit 9 —> 0. For a whatever 
small field h ^ this term goes to zero with 9 — > 0, and x_l — > 1/2. The 
transition to the regime, where the magnetic field exerts the influence on the 
susceptibility of a 2-dimensional antiferromagnet, is sharp due to the strong 
exponential temperature dependence of G ( 2. 22] ) and occures at the temperature 



2 ln(4/ft) 



(4.10) 



Note that the value of 9* is for h <C 1 much larger than the corresponding 
characteristic temperatures in the 1-dimensional case. The longitudinal suscep- 
tibility of the s.q. classical antiferromagnet has the form 



Xz 



1 



1 

D 



In 

7T 



l-G + h 2 /2J 7rl-G + h 2 /2 



+ 9 + 0(9 2 ) 



(4.11) 

where the addditional in comparison to (4.9) term is not very essential at low 
fields in contrast to the 1-dimensional case (see (3.4). The temperature depen- 
dences of x± and \z m the magnetic field obtained by the numerical solution 



of the equations (2.19) and (2.20) and the numerical integration in (4.4) are 



represented in figs. 6 and 7. 

The 1/D correction to the magnetization fl4.8| ) increases with approaching 
the magnetizatio n sat uration point h = 2 since here G = — 1 and P( G) ^> 1. 
But the formula (4.8) bec omes insufficient in this region, because in (4.4) the 
integral with ip' q (see (|4.3| )) becomes for G = — 1 comparable with the one with 
P'(G) due to the great longwavelength contribution, and one should use the 
complicated analytical expression for i p' f or x ~ 1 — G [01. In contrast, the 



quantity aipq in the denominators in (4.4) and (4.5) can be neglected in the 
whole field region, since in the low-temperature range aip = 9P(G) = 1— ttIq. 
In the field region above the saturation point (h > 2) at low temperatures 
m = 1 and in (jO)) G q = -G. Neglecting the terms BP' < 1 in AG flO] ), 
with the use of (2.23) for the total magnetization m = mp -I- mi/ D one gets 



i--^)GP(f; ): 



G = 



1 



h-1 



(/i>2,6»<1) (4.12) 



This result is the exact expression for the leading correction to the magneti- 
zation of a classical antiferromagnet in the spin-flip phase (H > 2|Jo|) in the 
low-temperature limit, which can be obtained independently with the help of 
the lowest-order spin wave theory. In the framework of the diagram technique 
for classical spin systems used here this corresponds to taking into account only 
the simplest diagram for the magnetization m with one integration over the Bril- 
louin zone (i.e. the one analogous to the 2-nd diagram in fig.l,b). The derivation 
of the formula (4.12) is trivial, because the ground state of the system has no 
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Fig. 6: Temperature dependences of the transverse susceptibility of the s.l. classical 
Heisenberg antiferromagnet for different magnetic fields in the 1-st order in 1/D. 
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spontaneous symmetry breaking and the magnon spectrum has a gap. However, 
with the approach to h = 2 in (4.12) G — > 1 and for low-dimensional systems the 
temperature correction to m diverges. In the region h < 2 the situation becomes 
complicated, and to obtain finite results for the thermodynamic quantities one 
has to take into account some infinite series of diagrams, wh ich is exem plified by 
the 1/D expansion described above. With the help of (4.8) and ( 4.12 ) one can 
write down the expressions for the magnetization of the 2-dimensional antifcr- 
romagnet on both sides of the magnetization saturation point h = 2 excluding 
a small intermediate region: 



In- 



2 ttD 2-h' 

2?r V D h-2 



Bin- -C2-/K 1 
a 

01n-</i-2<l 



(4.13) 



The latter results are analogous to those for the linear chain model ( |3.7|) that 
could be also obtained in the same way like here. The normalized susceptibility 
Xz of the 2-dimcnsional antiferromagnet is greater then 1/2 below the saturation 
point h = 2, too, and the maximal value of Xz is greater than that for the linear 
chain (see fig. 8). The latter can be explained by the fact that for a square lattice 
there is no competing negative contribution to m of the zeroth order in 1/D, as 
is the case for the linear chain (see (3.7)). 



5. Discussion 

In this article the 1/D expansion of the physical quantities of low-dimensional 
classical Z?-vector models in the whole range of temperatures and magnetic 
fields was developed, the results obtained being valid for both ferro- and an- 
tiferromagnets. For the calculation of the susceptibility and the field-induced 
magnetization of ferromagnets at low temperatures the method is, however, not 
very efficient, because these quantities are singular at T — > 0. In rcf. S was 
shown that at low temperatures the 1/D correction to the susceptibility of a 
2-dimensional ferromagnet becomes greater than its value in the zeroth order 
in 1/D, which means that D enters the argument of the exponentially great 



expression for x oc 1/(1 — G* ) (comp. ( 2.21 )). This is in accord with the results 
of the RG-approach of ref. |i7|j , but does not allow to obtain accurate results 
in the framework of the 1/D expansion. 

On the other hand, the 1/D expansion proves to be a very good instrument 
for the description of noncritical characteristics of low-dimensional magnets, 
such as the magnetization and susceptibility of antiferromagnets and the en- 
ergy and heat capacity of both ferro- and antiferromagnets. For the latter the 
zero- field results (identical in both cases) were obtained in ref. ]^], and their 
generalization to the case with the magnetic field with the use of the methods 
developed here makes no difficulties. The most remarcable feature of the results 



18 



3 _ 2j\ 


1 1 






1- 9=0.03 \ \ 




2- 9=0.1 \ \ 




3- 9=0.3 \\ 




4- 9=1 X 




e = T/7 MFA 




Square lattice classical 




Heisenberg antiferromagnet 



0.5 1 1.5 2 2.5 

h=H/\J \ 

Fig. 8: Field dependences of the longitudinal susceptibility of the s.l. classical Heisen- 
berg antiferromagnet for different temperatures. 



obtained with the help of the 1 / D expansion is that they describe the maximum 
in the temperature dependence of the zero-field antiferromagnetic susceptibility 
and its singular behavior at H, T — * 0. The former is the result of taking into 
account the diagrams with double integrations over the Brillouin zone, which 
was not made in any of the preceding theories. This means allowing for the 
wavevector dependence of the compact part of the spin-spin correlation func- 
tion (2.1) as well, or, in the other words, going beyond the Ornstein-Zernike 
form for Xk- It is worth noting, however, that at low temperatures (9 <C 1) 
the results simplify to some expression with only one integration over the Bril- 
louin zone (see (4.8)) which corresponds to taking into account the out-of-plane 
magnons having the field-dependent gap. 

An intriguing property of the l/D expansion is that it leads to the ex- 
act results for the noncritical characteristics of low-dimensional magnets at low 
temperatures. In particular, the physically expected limiting value of the an- 
tiferromagnetic susceptibility at T — > and H — is exactly recovered and 
the singularity of x{H> T) at H, T — > is entirely explained in the 1-st order 
in l/D. The latter it is unlikely to be essentially changed in the next orders 
in l/D. All the examples considered up to now suggest that the coefficients in 
the expansions of the noncritical quantities in powers of 9 are polinoms in 1 / D 
(sec rcf. Q). If it is true, then the up to now not available low-temperature 
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expansions of these quantities can be obtained with the help of the 1 / D expan- 
sion! Further it should imply that there is some method of derivation of these 
low-temperature expansions without using the l/D expansion. The search for 
such a method is planned for the nearest future. 

It would be very interesting to compare the results of the l/D expansion 
with results obtained by other methods. In particular, for the energy of a 
square-lattice classical Heisenberg magnet the MC simulations were made by 
Shenker and Tobochnik (see the comparison in ref. Q), but the antiferro- 
magnetic susceptibility was simulated by various researches only for a quantum 
model with S = 1/2. As concerns the 2-dimensional model with D = 2, the 
l/D expansion cannot, naturally, describe the Kosterlitz-Thouless transition, 
which occurs in this system. But one can expect that the general features of 
the temperature dependence of the antiferromagnetic susceptibility in magnetic 
fields described by the l/D expansion are inherent in this model, too. More- 
over, in the magnetic field the behavior of the antiferromagnetic model with 
D = 2 should simplify, because the magnetic field lifts the spontaneous sym- 
metry breaking and induces the gap of spin fluctuations. In this case at low 
temperatures it is enough to take into account only the lowest order diagram 
of the spin- wave theory that is naturally contained in the l/D expansion in 
the 1-st order in l/D (comp. (4.12)). It should be stressed that for the model 
with D = 2 the effects in the temperature and field dependence of the antifer- 
romagnetic susceptibility discussed in this paper should be the most strongly 
revealed. For the Heisenberg (D = 3) antiferromagnet the behavior in the mag- 
netic field can be more complicated, than for D = 2. As we have seen, at low 
temperatures even small magnetic field forces the spins to lie perpendicular to 
it. This decreases the effective number of spin components from D = 3 to D = 2 
and should lead to the Kosterlitz-Thouless transition with disappearance of the 
gap. But it should not change essentially the results for the susceptibility in 
this region, since in the expression for the 1 /D correction to the magnetization 
(1.5 ) enters the gap of the out-of-plane spin waves, which cannot disappear due 
to the Kosterlitz-Thouless transition. 

It should be also mentioned that the D-component vector model || consid- 
ered in this article can be generalized for hamiltonians with the spin anisotropy. 
For example, one can consider the so-called n — D model [flOf| , in which only n 
from the total D spin components are coupled by the exchange interaction. In 
this sence the x — y model (D = 3, n — 2) is something different from the plane 
rotator model (D = n = 2). It should be noted that the quantitis n and D 
play the different roles: the well-known expansion of the critical indices of 3- 
dimensional systems is the expansion in 1/n, and the expansion developed here 
for low-dimensional systems is the l/D expansion. The results of the present 
article can be generalized for the n—D model, as well as for more general models 
with anisotropic spin interactions. 
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Appendix 



The 1/D diagrams 

The additional diagrams constituting Aaa (k) in the expression f or th e com- 
pact part of the transverse spin-spin correlation function A act (k) (2^) repre- 
sented in fig. 2 have the following analytical form: 



A(2) 

AO) 

aa 

AW 

aa 

AS(k) 

AL 6 2(k) 

AL 7 Q +r) (k) 

AL 8 2(k) 



A - A 1 f dq R~ 

^aazz^z — ^aazz~^ ^j~^2, 

- A act/ 3 i a 77 A / 3 /377 Vojj^^VqVq 



A 2 



(27T) 

*,3/377^/3/3zA 77 

rfq 



7:A QQ/ 3 / 3 77 A^ / 3 z A 77 ,2 va( T ^ r ^f3J Z qV^ 



(2tt) 

(2tt)^ 



(2tt) c 



r/q 



c/q 



(A.l) 



A.... • .'A,,,,- -.A ,'i.A- , f o/ ^ d /3 Jk-q/3 J zq V^ 



Here A QQ( g^ = A / 3^ 77 , A aQ z = Ag^, etc., are the renormalize d mu lti-spin 
cumulants with a ^ /3 =/= 7 (no summation over /3 and 7 in (A.l)) given 
by the formulas analogous to As the diagrams (A^l) should be calcu- 

lated only in the 1-st nonvanishing order in 1/D, one can use for the renor- 
malized transverse interaction line /3J q (see fig.l,c) the simplified expression 
/3J q = /3J q /(l — A aa (3Jq), where A aa = {9/D)G and G corresponds to the 
spherical model (see ( [2.20 )). The renormalized longitudinal interaction /3J zq 
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(see fig.2,b) is given by 



BJq 



i-(A K uLm 



(A.2) 



and the renormalized 4-spin correlation line Vq (see fig.2,c) reads Vq = Vq/(1 — 
haappVq), where 

V« = jV O J^f3J p 0j^ p (A.3) 

is the unrenormalized 4-spin correlation line, the factor D (or D — l, which plays 
no role here) in (A.3) resulting from the summation over the spin-component 
indic es 3 and 7 in the diagrams. Calculating now the renormalized cumulants 
A in ( A.l ) in the lowest order in 1/D by the pass method, one gets (|] 



A Qa /3/3 



DJ (2y) 2 (2y-l)' 
with y = 1/ (9G) and, additionaly, 



A 



A 2Z 



A r 



A QQ /3^^ 2 



aQ/3/377 



A, 



2(3y-l) 
D; (2y)3(2y-l)3 

(A.4) 



(A.5) 



where £ = (D /9)(h + vm^) = Dym^. With the use of these results for the corre 

(D 2 /e)A£l D) (k) 



spending contributions into A^ ' ; = limn 



A (D 

A (2) 

A< 3 > 

A ( 4 ) 

AS(k) 
A( 6 + 7+7 ' +8) (k) 



Sl/D) = 

vG 2 
"(2tf-l) 

-2G 5 3y_1 



1 ^2 % - 1 

1 - 4m y 



(2i/ lY V °J('2m' 
r,;/w (2y-l)4^(27r)^ 



= —%vm1yG 



{2y-lf 
dq 

¥>q^q 

dq - 



l>0 



, one obtains 
dq ~ 



(27T) 



' I'O 



A 2q <^q 



wG (2^ry 



(A.6) 



2zyG 4 



(2y-l)5 
dq - 



2tt) 



(2y- 1)^(27^ 

4mlyG 3 f dq 
(2y-l) 2 ^o/(^I^- q A 2qLq 

where <^ q = <^ q L q = <^ q /(l + aipq), a = G 2 /(2y - 1), 



k-q^q 



dp ~ ~ 



(27T) 



d ApA p _ q ; 



A q = J q / Jo = 



An 



1 - vG\ c 



(A.7) 
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The expression for the renormalized z-interaction line A zq = J zq /Jo can be 
written in the form 



Further simplifications leading to the results listed at the end of the Section 
2 can be archieved if one expresses (fq through r q ( 2.32| ) with the use of the 
relation (2y — 1)(1 + a<Pq) = r q + 1m\y. 
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